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The JCMT Gould Belt Survey: SCUBA-2 observations of radiative 
feedback in NGC1333 
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ABSTRACT 

We present observations of NGC1333 from SCUBA-2 on JCMT, observed as a JCMT 
Gould Belt Survey pilot project during the shared risk campaign when the first of four arrays 
was installed at each of 450 and 850 microns. Temperature maps are derived from 450|^m 
and 850|i.m ratios under the assumption of constant dust opacity spectral index [3 = 1.8. 
Temperatures indicate that the dust in the northern (IRAS 6/8) region of NGC1333 is hot, 20- 
40 K, due to heating by the B star SVS3, other young stars in the IR/optically visible cluster, 
and embedded protostars. Other luminous protostars are also identified by temperature rises 
at the 17" resolution of the ratio maps (0.02 pc assuming a distance of 250 pc for Perseus). 
The extensive heating raises the possibility that the radiative feedback may lead to increased 
masses for the next generation of stars. 
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1 INTRODUCTION 

Gas temperature plays a role in c ontrolling st ar formation in dense 
cores through thermal support ( |jeanslll902h . The thermal Jeans 
length may ultimately control the sub-fragmentation of cloud cores 
once turbulent motions have decayed and cores become coher- 
ent (Goodman et al. 1998; Myers 1983). Heating of protostellar 
cores suppresses sub-fragmentation (i Bate 2009; Offner et al. 2009; 
lUrban. Martel & Evansll201(]| ; IChabrier & Hennebelle..201 1.) and is 
potentially an important factor in the formation of massive stars 
(along with magneti c fields, co mpetitive accretion, and other fac- 
tors; see Hennebell e & Commer con 2012 for a review). 

Observationally, several techniques exist for measuring tem- 
perat ures of molecular gas, inc l uding C O excitation temperatures 
(e.g. ICurtis. Richer & Buckl3 |2010|; ITadd. Mvers & GoodmanI 
1 19941) and the ammonia ladder (e.g. Huttemeister et al. 1993). The 
dust colour temperature Td provides an alternative estimate of the 
gas temperature: at high molecular hydrogen densities > 
IC*'^ cm~^, where gas-grain collisions dominate the gas heating, 
the dust and gas temperatures are well coupled tGoldsmith 200 L) . 
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Dust colour temperatures rely on fitting dust continuum fluxes to 
an opacity-modified blackbody spectrum, a method which is being 
applied very successfully to make temperature and column density 
maps of nearby star-forming regions from Her schel data based on 
fits from 160 ^m to 500 |i.m jAndre et alj201(]t) . 

In this paper we apply dust temperature methods to 450p.m 
and 850 p.m submillimetre (submm) dat a from SCUBA-2 on the 
James Clerk Maxwell Telescope (JCMT; [Holland et aTll2006h. us- 
ing a similar technique to that developed for SCUBA jKramer et alj 
I2OO3I) . There are benefits in using long-wavelength mm/submm 
data for dust temperatures (as long as the data points are not en- 
tirely on the Rayleigh- Jeans tail). When temperature varies along 
the line-of-sight, dust continuum fluxes are emission-averaged and 
hence weighted towards the hottest dust. Small fractions of hot dust 
can dominate the mid-infrared (mid-IR) even though they contain 
only a small fraction of the total dust mass , as is evident from pro- 
tostellar spectral energy distributions (e.g. lEnoch et alj|2009l) . For 
sources with a simple geometry, this temperature variation issue 
can be circumvented by radiative transfer modelling, but this is not 
a realistic option for complex star-forming regions. Also, in high 
column density regions, the dust emission at mid-IR wavelengths 
can become optically thick, weighting the results towards the front 
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side of the cloud and low-opacity photon escape routes. The 14" 
beam full width half maximum (FWHM) of JCMT at 850 ^m 
provides an improvement in angular resolution over Herschel at 
500 [Lva (37" FWHM), which is advantageous in dense clusters 
and when determining the relative location of heating sources and 
cores. 

Here we apply the dust colour temperature method to 
the young cluster NGC1333 in Perseus. NGC1333 is a young 
cluster well studied in the mm/submm ('Sandell & Rnee^ '200 1'; 
Hatchell et al. 2007; Kirk, Johnstone & Tafalla 2007; Enoch et al. 
20061 ; ICurtis, R icher & Buckle 2010; Curtis etalj |2010l) ; for gen- 
eral background see .Walawender et al.. (2008.) . 



2 METHOD 

2.1 Observations and data reduction 

NGC1333 was observed on February 21-2 2 2010 and March 4th 
and 7 th 2010 as a JCMT Gould Belt Survey jWard-Thompson et alj 
l2007h pilot project in the SCUBA-2 shared risk (S2SRO) cam- 
paign, during which only one array out of the final four wa s 
available at each of 850^m and 450 vim ( [Holland et"al] l2006h . 
These observations were among the first made with SCUBA-2. 
The weather conditions were good, JCMT Band 2 with 225 GHz 
opacity T225 = 0.05-0.07. Two overla pping 15' diameter circu- 
lar regions ('PONG900' mapping mode; |jenness et ^1201 ih were 
targetted, with map centres offset 10' to the north and south from 
RA(2000) = OS*" 29'" 00:0, Dec. (2000) = 31° 18' 00" . These 
fully-sampled 15'-diameter regions were surrounded by 4' borders 
that were scanned at lower sensitivity by decreasing numbers of the 
array pixels. 

Molecular lines contribute to the signal in the SCUBA-2 
850 [im continuum bandwith, with the most signifi cant contam- 
inatio n at 345 GHz due to the CO J = 3-2 line jPrabek et al] 
The CO contamination was removed during the 850 ^tm 
reduction usin g the HARP ^^CO 3-2 integrated intensity map 
jCurtis, Richer & Buckle 2010) converted to a CO contamination 
map using a contamination factor of 0.66 mJy/beam/(K km s~^) 
appropriate to Ban d 2 weather with an average T225 = 0.06 
jPrabek et al.l2012t) . The 850 jxm map was re-reduced in June 20 1 2 
using th e development ver sion of the data reduction package 
SMURF jjenness et al.l201ll. version of 2nd May 2012) with astro- 
nomical signal (AST) masking based on earlier versions of the re- 
ductions and gridding to 6" pixels. The CO contamination map was 
injected during the 850 i^m reduction process as a fake source with 
a negative multiplying factor of — 1/FCF where FCF is the flux 
conversion factor (see below). In this way, the CO emission was 
subtracted from the timeseries fluxes during the reduction, match- 
ing the spatial filtering inherent in SCUBA-2 map reconstruction 
and resulting in an 850 \im map which is corrected for CO line 
contamination. The 450 (xm map was re-reduced in May 2012 us- 
ing the same SMURF version, with AST masking, without CO con- 
tamination removal, and gridded to 4" pixels. Both the 850 \im and 
450 [im maps are insensitive to emission on scales greater than the 
on-sky footprint of the array (4' for the single-array S2SR0) due 
to common-mode subtraction of atmospheric emission across the 
array. 

Maps were calibrated by multiplying by flux conversion 
factors (FCFs) of 653 ± 49 Jy/beam/picoWatt and 511 ± 31 
Jy/beam/pW at 850 |J.m and 450 \im respectively, d erived from the 
six observations of the submm calibrator CRL618 jPempsev et alj 



l2012h . The average beam PWHM were measured to be 14.4" ± 
0.15" (850 |J.m) and 9.8" ± 0.23" (450 ^.m) consistent with the 
observatory standards of 14.2"/9.4". 

Of the 18 maps observed, 12 were selected for mosaicking, 
six each for the north and south centres, with six rejected due 
to telescope tracking errors. The resulting mosaics have pixel-to- 
pixel 1(7 noise levels of 7.8 and 130 mJy per 14.4"/9.8" beam 
on 6"/4" pixels at 850 (i.m/450 [im respectively (an average 
over the north and south regions, which have noise levels -1-4% 
higher and lower respectively), roughly a factor two deeper than 
SCUBA (Sande ll & Kne3l200ll ; [Hatchell et al.ll2005l) . Calculating 
the equivalent point- source sensitivity at a resolution equal to the 
beam area, these observations meet the JCMT Gould Belt Survey 
3 mJy target at 850 \im while achieving 50 mJy point source sensi- 
tivity at 450 i^m. The central, overlap region of the map has a still 
lower noise level by roughly \/2. The pixel-to-pixel noise estimates 
do not completely represent the noise levels in larger-scale fluctua- 
tions which were measured from histograms of the pixel values to 
be 14 mjy/beam and 150 mjy/beam respectively. 



2.2 Ratio and temperature maps 

For optically thin emission, the flux ratio between the two wave- 
lengths 450 [xm and 850 |J.m depends solely on the dust opacity 
and temperature. The dust opacity is stan dardly para meterised at 
long wavelengths as a power law k oc ly^ ( [Hildebrmdi 1 983.V The 
450/850 flux ratio is then 

S450 ^ / 850 V"^^ / exp(/tc/A85ofcrd) - l \ 
Ssso ^ \^50 J \exp{hc/\45okTA) - 1 J 

where Td is the dust temperature and Agso and A450 are 850 [im 
and 450 \im wavelengths respectively. 

Higher 450/850 flux ratios can be produced either by higher 
dust temperatures, a higher opacity index /3, or both. Separating 
the two requires data at additional wavelengths, e.g. a joint analysis 
with the Herschel 250 [im and 500 vim datasets, which is planned 
for future work beyond the scope of this Letter. Higher dust tem- 
peratures are caused by heating, whereas grain growth can change 
P by raising the typical sizes and absorption wavelengths of grains. 
P typically lies in the range 1 to 2, depending on the type and 
dumpiness of grains. Increases in B hav e been observed in pro- 
toplanetary disks ( e.g. [Ubach et alJIioT^ and low-density clouds 
dMartin et alj20I j) but to date there is little knowledge of changes 
in pre/protostellar cores, a situation that should improve rapidly 
with availability of SCUBA-2 and Herschel data. There is a sugges- 
tion that icy mantles increase towards the inner parts of the Class 
sourc e B335 though not enough to significantly affect /3 dPotv et al] 
l20Tdh . 

The flux ratio map was created from the 450 [im and CO- 
subtracted 850 vim maps, which were scaled to Jy/pixel with any 
negative fluxes removed. The 450 vim map was convolved with a 
2-dimensional Gaussian representing the 850 vim 14.4" primary 
beam, and the 850 vim map similarly convolved with a 450 vim 9.8" 
primary beam to match angular resolution. The cross-convolution 
procedure produced more consistent matched beams than simply 
convolving the 450 vim map to the 850 vim resolution and was 
needed for the S2SR0 data to avoid artefacts in the ratio map. The 
flux ratio map was created by dividing the resolution-matched, re- 
gridded 450 vim and 850 vim maps with values in each map below 
5(T blanked. For a fixed value of the dust opacity /?, Equation [T] 
was used to convert the ratio map into a temperature map, with an 
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additional cut of temperatures with variance > 30 K applied. The 
effective angular resolution of the ratio and temperature maps is 
17.4". The systematic uncertainty in the ratio maps due to calibra- 
tion at each wavelength is 15%. The large-scale filtering inherent 
in S2SR0 maps adds edge effects and systematic offsets to the ra- 
tios which will be quantified by future modelling. These effects are 
smallest where bright compact structures dominate the fluxes, as is 
the case for most of NGC1333. 



3 RESULTS 

The S2SR0 850 ^m map of NGC1333 is presented in Fig.Ea) 
as a finding chart. SCUBA cores a re lab elled with their HRF2005 
numbering ( iHatchell etal] l2005l l2007l hereafter HRFnn) and 
rotostars identified from Spitzer mid-IR de tections are marked 
Gutermuth et al. I I2OO8I : |j0rgensenetai]|2OO6h . Further identifica- 
tions of sources in NGC1333 can be found in ISandell & Kne3 
( l200lh . Figure Hb) shows the 450 \mi map of this region. The 
450 \Lm dataset from S2SR0 is sensitive and stable enough for the 
filamentary structure to be seen clearly, as it is not dominated by 
the large-scale noise artefacts which limited analysis of the SCUBA 
450 ^m data ( Sandell & Knee 2001; Hatchell et al. 2005). 

In Fig. [TJc) we plot the dust temperature in NGC1333 as- 
suming a constant dust opacity spectral index of /? = 1.8. This 
choice of /3 is consistent with the p opular OH5 dust model (/3 — 
1.85; lOssenkopf & Henning|[l994l) , with recent results for dense 
cores from Planck and Herschel jjuvela et alj I2OI ll ; IStutz et al] 
and with compar isons of SCUBA and near-infrared opac- 
ities dShirlev et alj|201lh . Although variation in the 450/850 ratio 
may also be due to variation in (3 (to be investigated in future work), 
as NGC1333 contains several luminous protostars and a reflection 
nebula, it is reasonable as a first approximation to assume that the 
variation is dominated by temperature. A higher value for /3 would 
result in lower temperatures for the same ratios: temperatures of 
10 K, 20 K and 40 K at /3 = 1.8 fall to 8.5 K, 14 K and 20 K at 
13 = 2.2. 

Filaments and ambient cloud material generally show dust 
temperatures of 10 K (corresponding to 450/850 ratios of around 
4; ). These temperatures are typical of the cold dense interstel - 
lar medium jSvans et 'Zl l200ll ; lshirievy Evans & Rawlingsll2002h . 
Towards the protostars, beam-averaged temperatures rise to 12- 
15 K, consistent with SCUBA-based models of low-mass Class 
and 1 protostars dShirlev, Evans & RawUngsl I2OO2I ; | Young et al 



20031) and with the Herschel analysis of CB244 Istutz et al 



2010). Some of the more luminous protostars show notably higher 
dust temperatures as discussed below. The temperature does not, 
however, rise towards all protostars identi fied by Spitzer mid- 
IR detections or outflows 1 J0rgensen et al.l 2006; Hatchell et al. 



20071; iGu termuth et al.' ' 20081 ; Inatchell. Fuller & Richer! i2007l ; 



Hatchell & Dunham 2009). Low-luminosity protostars have simi- 



lar temperatures of 10-15 K to starless cores. The falling 450/850 
ratios towards the map edge are almost certainly an artefact of the 
filtering of the large-scale background in the S2SRO map recon- 
struction, and not a physical effect due to falling /3 or tempera- 
ture. Edge reconstruction should improve with data from the full 
SCUBA-2 array, which is sensitive to spatial scales a factor of two 
larger. 

The most obvious feature in the temperature map is the warm 
northern region containing cores HRF45, 47, 54, 56, and 66 that 
shows widespread dust temperatures above 20 K. This region lies 
to the south of the optical cluster, and its dust is presumably 



heated by the stars in the optical cluster as well as the embed- 
ded protostellar population. The opt ically revealed cluster stars 
are dominated by th e B5 star SVS3 l lStrom. Vrba&Stron]| 19761; 
IStraizvs et all l2002t) which lies on the edge of the dust to the 
northeast of HRF54, as marked in Fig. [Ha). SVS3 is a 138 Lr 
binary with an F2 companion dCormellev. Reipurth & Tokunagal 
.200 81) ; the optically-visible NGC13 33 reflection nebula and as- 
sociated infrared source (IRAS 8; Ijeimings et al.l 1 19871) is ev- 
idence its proximity to the dust cloud. The other main heat- 
ing sources are embedded p rotostars. Core HR F45 is heated by 
SSTc2d J032901. 6+3 12021 ( lEvans etalj 12009) as suppo rted by 
its hi gh temperature and identification as IRAS 6 (Jennings et alj 
Il987h . Protostar SSTc2d J032907.8+312 157 in HRF56 is also 
classified as a 27 Lp) luminosity source jGutermuth et al.l I2OO8I ; 
lEvans et alj|2009l) . This luminosity is likely to be partially due to 
accretion and partially reprocessed flux from SVS3. The dust at 
this position is warm (~ 20 K), though not as warm as neigh- 
bouring HRF54 heated by SVS3 to form IRAS 8. The ridge 
containing HRF66 has elevated temperatures; a protostar embed- 
ded in HRF66 could explain some of this heating but a second 
source would also be necessary to heat the ridge to the south- 
west. Sources HRF47 and 53 are also warm (~ 20 K) but do 
not show temperature peaks. Core HRF4 7 contains a Spitzer- 
identified p rotostar, J032859.56+312146.7 /Jer gensen et aUboOd ; 
lEvans et "alj 2009), but the closest Spitzer detection to HRF53 is 
SSTc2d J0329 05 . 1 8-)-3 1 2036.9 to the southeast where th e dust tem- 
perature rises I Gutermuth et alJl2008l ; lEvans et alj20"09l) . 

An insight into the location of the heating sources can be 
made by comparison with ^^CO excitation temperatures Tex (CO) 
dCurtis Richer & Bucklj|2010l. Fig. 13). The ^^CO traces warm 
gas on the front side of the cloud (as the ^^CO is generally op- 
tically thick) but the submm dust emission is optically thin and 
traces heating through the cloud. Towards SVS3, the ^^CO tem- 
perature peaks sharply at over 40 K indicating heating of the front 
of the cloud. The dust temperature rise is more extensive, suggest- 
ing additional heating from embedded sources or sources behind 
the cloud, including SSTc2d J032905. 18+3 12036.9. Towards the 
HRF45 and 47 cores, the dust and ^^CO excitation temperatures 
are similar (~ 40-45 K), indicating that the heating sources lie 
near the front side of the cloud. 

Further to the north, another B star, BP -1-30 54 9 or 
NGC1333 IRAS 9 classified as B8 llenning s et al] Il987l) . lies 
close to HRF57 and 67 in projection which may explain raised 
dust temperatures of 15-17 K in these northern filaments. A sig- 
nificantly higher CO temperature, Tex (CO) over 40 K, suggests 
that BD +30 549 lies in front of the cloud. Several protostars in 
other parts of the cloud are identifiable by their raised 450/850 ra- 
tios and temperatures. NGC 1333 IRAS 7 (HRF46) shows raised 
Tex (CO) ~ 35 K and raised Tdust ~ 20 K to the northeast of the 
dense core. This could be due to the embedded source lying towards 
the front side of the cloud, or more likely in this case, h ot CO lining 
a pro tostellar outflow cavity opening to the northeast jCurtis et alj 
[20T0I) . 

Deeply embedded sources or sources behind the cloud show 
the opposite effect with Td > Tex (CO). These include the 
SVS13 ridge (HRF43) and NGC 1333 IRAS 2 (HRF44). The 
HRF49 condensation shows a particularly strong dust temper- 
ature gradient from lOK in the southwest to over 40K at the 
northeast head of the conden sation, heated by infrared source 
SSTc2d J032837. 1+311331 ( IGutermuth e t al. 2008; Enoch et alJ 
1 2009) . ^^CO excitation temperature maps fCur tis, Richer & Bucldg 
l20I(i Fig. 13) show no evidence for raised gas temperatures in this 
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Figure 1. Left to right: (a) CO-subtracted S2SR0 850 ma p of NGC1333. Contours are at 0.02,0.1,0.5, 2.5 Jy/beam. The bl ack Une marks t he CO subtraction 
region. Labels identify SCUBA sources ( Hatchell et al.'2005') and the two B stars. Markers identify embedded protostars from lHatchell et al] i2007 l) (red circles 
for Class 0, orange circles for Class I), J0rgensen et al. (2006) (yellow +) andfOutermuth et al. |200§J (yellow x). (b) S2SR0 450 \im map in Jy/beam. (c) 
Dust temperature map at /3 = 1.8 in units of Kelvin. Contours are at T^j = 12, 20, and 30 K. 



region indicating that the front of the core from which the optically- 
thick ^■^CO escapes is cold and the star and the hot dust lie behind 
the dense core. 



4 DISCUSSION 

Our observations show that star formation has raised 450/850 ratios 
in NGC1333. The absolute values of temperature in Fig.[T]depend 
on our assumption of /3 = 1.8, and increases in /3 may account for 
some of the 450/850 ratio rises towards dense cores. However, an 
explanation based on opacity alone would require extreme values of 
/3 ~ 3 towards the reflection nebula and luminous SVS13, whereas 
dust heating provides a simple explanation for ratio rises in these 
regions. 

According to theory, heating increases thermal pressure which 
acts t o stabilise cloud material against further fragmentation dJeansi 
fl903) . The thermal Jeans length and Jeans mass scale as \/T and 
T^/^ respectively. Cores may be heated internally, potentially sup- 
pressing fragmentation resulting in more massive stars, or exter- 
nally by earl ier generations of sta rs, as recently investigated in 
G8.68-0.37 ( lLongmoreetalj201lh . In NGC1333, there is no evi- 
dence that the northern region has fragmented with larger core sep- 
arations than the rest of the region to date, presumably because the 
l'-1.5' core separation sc ales have been set by turbule nt, rather 
than thermal, support (e.g. IChabrier & HennebelldfloTTb . The lo- 
cal temperature rise in NGC1333 N due to the existing optical 
cluster and embedded protostars may, however, now inhibit the for- 
mation through fragmentation of further gravitational cores within 
the northern region. At the average molecular hydrogen density of 
NGC1333 N, riH, — 4 x 10* cm~'^, the Jeans mass increases from 



0.5 to 4 Mq as temperature rises from 10 to 40 K. If no further 
cores form, the existing protostars will accrete the remaining grav- 
itationally bound gas (minus a fraction removed by outflows). The 
existing cores will produce more massive stars than they would 
have done in the absence of external heating. Hence the radiative 
feedback from the existing cluster may influence the development 
of the next generation of stars. 

There is possible evidence in the data here of the effect 
of external heating in promoting massive star formation in core 
HRF56. The embedded protostar SSTc2d J032907. 8-1-312157 lies 
only 0.05 pc in projection from the B5 star SVS3 and on the 
edge of the optic al nebula. If its current luminosity of 27L0 
tevans et alll2()09h is internally generated (rather than reprocessed 
emission from SVS3), and with a few M© envelope still to ac- 
crete, HRF56 could be forming a second B star. Its lack of frag- 
mentation could be influenced by heating from nearby SVS3. The 
nearby HR F54 core may be affect ed similarly in the future. We 
agree with iLongmore et al.l ( 1201 1[) that this is a somewhat un- 
satisfactory circular argument as it requires a pre-existing mas- 
sive star (SVS3) to provide the feedback, rather than a cluster of 
low-mass stars. Further to the south, HRF45/IRAS 6 has a cur- 
rent luminosity of 15L0 and a few Mq envelope so could also 
reach B-star status. This source, however, does not lie close to any 
such luminous star and a double radio detection indicates that it 
has already fragmented ( Rodriguez, Anglada & Curiel 1999). Any 
further stars formed are likely to be typical AFGK-type clus- 
ter members since the remaining cores are low mass, at most a 
few Mq each. (Note that the total mass in NGC1333 N is only 
~ 2OM0 , assuming an average dust temperature of 20 K. Previous 
core mass calculations assumed fixed t e mperatures ( Hatchell et alj 
I2OO7I : iKirk. Johnstone & Tafallal |2007| : ISandell & Kned l200lh or 
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temperatures from IRAS jSandell & Kneel200ll) . A future step will 
be to recalculate the masses of the NGC1333 N cores based on the 
dust temperatures and fluxes from the full SCUBA-2 Gould Belt 
survey data, which will have better spatial response and stabler cal- 
ibration than the S2SR0 data presented here). 

Finally, in HRF54/IRAS8, radiative feedback by SVS3 may be 
promoting formation of an early B-type star with potentially dev- 
astating consequences for NGC1333 through widespread dissoci- 
ation of the molecular gas and a halt to further star formation ac- 
tivity in the region. Possibly, the high-mass stars currently missing 
from the IMF of Pe rseus and other young low-mass star-forming 
regions HE vans et alj[2009.) only form once radiative feedback from 
the existing clusters reaches a sufficiently high level. A further con- 
sequence is that massive stars are likely to form relatively late in 
molecular clouds, in the centre of clusters where the gas heating is 
greatest. 



5 CONCLUSIONS 

SCUBA-2 S2SR0 observations show significantly raised dust tem- 
peratures around several luminous protostars in NGC1333, includ- 
ing SVS13 and IRAS 2. Assuming a dust opacity spectral index 
/3 = 1.8, the region on the southern edge of the optical reflec- 
tion nebula is extensively heated to temperatures > 20 K, reaching 
40 K close to the B star SVS3 and the embedded protostar IRS 6. 
A significant fraction of the energy input is from the optical cluster, 
and specifically the B star SVS3 illuminating the reflection nebula. 
Heating by the existing cluster is potentially biassing the star for- 
mation in the north of NGC1333 in favour of more massive stars. 

NGC1333 will be observed again with the fully- 
commissioned SCUBA-2 as part of the JCMT Gould Belt 
Survey. 
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